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u i ABSTRACT 

<: 

CT) ' Recent observations with the Hinode Solar Optical Telescope display an active 

, region prominence whose fine threads oscillate in the vertical direction as they 

move along a path parallel to the photosphere. A seismological analysis of this 
event is carried out by taking advantage of the small radius of these structures 
compared to the total length of magnetic field lines, i.e. by using the thin tube 
approximation. This analysis reveals that the oscillatory period is only slightly 
modified by the existence of the flow and that the difference between the period 
CN ■ of a flowing thread and a static one is below the error bars of these observations. 

■ Moreover, although it is not possible to obtain values of the physical parameters, 

a lower bound for the Alfven speed (ranging between 120 km s _1 and 350 km s _1 ) 
{Nj ■ is obtained for each of the threads. Such Alfven speeds agree with the intense 

CO ■ magnetic fields and large densities usually found in active region prominences. 

o : 

OO . 

Subject headings: Sun: prominences — Sun: oscillations — MHD — waves 

> : 

d | 1. Introduction 

Solar prominences are cold (T ~ 10 4 K) plasma structures embedded in the much hotter 
solar corona. The prominence material is much denser than its surroundings and, since it 
stays suspended in the corona longer than the free-fall time, some force must counteract the 
downward force of gravity. Such an action is thought to be of magnetic origin and one of the 
prevailing ideas is that the prominence plasma lies on a dip of the magnetic field lines and that 
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the field line curvature provides with the required upward force (e.g. lAnzer fc Heinzell 120071 ). 
A common feature of many prominences is that when observed with sub-arcsec resolution a 
complicated internal structure manifests: prominences are composed of many thin, parallel 
threads embedded in material which does not emit or absorb in prominence spectral lines. 
The thread plasma often fl ows along the m agnetic field at speeds comparable to or larger 
than the local sound speed (ILin et al.ll2005l ). 



Prominence temperatures can usually be well determined by spectroscopic means, but 
other parameters are no t so well-known or change wide ly between different prominences. 
Tandberg-Hanssenl (119951 ) and iPatsourakos fc Viall (120021 ) give a compilation of various den- 
sity determinations in prominences and show that this parameter varies by at least two 
orders of magnitude (namely from 10 9 cm -3 to 10 11 cm -3 ). Magnetic fields show the same 
dispersion of values and range from a few G to 20-30 G i n quiescent prominences (e.g. 
Bommier et al.l 1 19941 ; iMerenda et al.l 120061 ; iGunar et al.l 120071 ) and even higher values in ac- 
tive region prominences. 

Given the difficulties in determining prominence parameters, prominence seismology 
has been put forward as an alternative metho d to probe the nature of these structures (see 
Oliver fc Ballesterl |2002| ; iBanerjee et al.l 120071 . for some reviews of this topic). The present 
work describes a seismological application to data obtained with the Hinode Solar Optical 
Telescope (SOT). 



Observations 



The Hinod e SOT observations analyzed in this work have been described in detail by 
Okamoto et al.l (120071 ). The data show an active region limb prominence composed of a 
myriad of thin horizontal threads that flow parallel to the photosphere. Some of these 
threads display vertical oscillations as they flow and these oscillations are synchronous, i.e. 
each thread moves entirely with a constant p hase. The main featur es of the moving threads 
and their oscillations are given in Table 1 of lOkamoto et al.l (120071 ) and summarized in our 
table [lj so we just remark that these are thin structures whose width and length range 
from 360 to 660 km and from 1700 to 16,000 km, respectively. The period and the velocity 
amplitude of the vertical oscillations are in the range 135-250 s and 8-22 km s _1 , while 
the horizontal flow velocity varies between 15 and 46 km s _1 . A prominence thread is a 
cold plasma condensation that occupies a segment of a much longer magnetic tube. In the 
present case it i s not easy to directly measure the length of such magnetic tubes, although 
Okamoto et al.l (120071 ) estimate that the wavelength of the oscillations is at least 250,000 km 
and so the minimum length of the magnetic tubes is 125,000 km. Another estimation of this 
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quantity can be obtained as follows: the oscillating threads lie parallel to the photosphere 
across the field of view of the SOT (which is 80,000 km wide), so that the full length of 
magnetic field lines is at least 80,000 km plus twice the vertical distance between the threads 
and the photosphere, which from table [T] is at least 12,400 km. All this amounts to a 
minimum magnetic tube length of 100,000 km. 



3. Theoretical analysis of thread oscillations 

Since every thread displays a distinct period of oscillation it is reasonable to assume that 
its dynamics can be studied by isolating each of these structures from the others. Moreover, 
curvature and the non-linear terms in the MHD equations are neglected and the low-/? limit 
is used. Thus, it will suffice to consider the transverse vibrations of a straight magnetic tube 
embedded in the corona and with a flowing segment consisting of dense, prominence plasma. 

We then adopt an equilibrium model (Figure [T]) consisting of a cylinder of length 2L 
embedded in a uniform corona with density p c . The plasma in the cylinder flows with speed 
Vo and consists of a dense, middle section representing the observable thread in the Hinode 
SOT data, and an evacuated part of coronal material. These two regions have uniform 
densities p v and p e . The plasma flow inside the magnetic tube results in the advection of the 
density structure with speed v . We assume that at t — the thread occupies the center of 
the magnetic tube, so the density along the tube can be expressed as 



p(z,t) 



p p if \z — v t\ < W, 
p e if \z - v t\ > W, 



where 2W is the length of the thread and the z-axis points along the tube, with z = the 
position of its middle point and z = ±L its end points. The magnetic field, Bo, is the same 
inside and outside the cylinder and parallel to its axis. 

The dynamics of the threads is studied taking into account that the ratio of the thread 
radius to the total length of field lines is a very small quantity, below 10~ 2 for all threads, 
so the thin tube approximation can be considered. A simple way to derive the equation 



govern ing the transverse oscillations is to follow the procedure of iDymova &: Ruderman 



(120051 ). who studied the normal modes of a structured magnetic cylinder, like the present 
one, for vq = 0. It is only necessary to include the two ingredients that are absent in 
that study: the flow inside the magnetic tube (which leads to the substitution of d/dt by 
d/dt + v • V = d/dt + vod/dz inside the cylinder) and the fully time-dependent situation 
(which prevents us from considering the time dependence exp(ioot)). It is worth to mention 
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that the terms coming; from the equilibrium flow can be ignored because, such as noted by 



Dymova Sz Rudermanl (120051 ). inside the cylinder the terms with derivatives along the tube, 
d/dz, are much smaller than those with radial or azimuthal derivatives. We then arrive to 
the remarkable result that the flow does not appear explicitly in the equations although it 
is present through the density; see equation JT]). Finally, a single differential equation that 
governs the transverse vibrations of the cylinder in the thin tube limit is obtained, 

d 2 u 9 . , d 2 u , , 

w -4 { ,,O w = o, (2) 

where u(z,t) is the transverse velocity component at the tube boundary, or in other words, 
u(z,t) is the velocity component responsible for the observed later al displacement of the 



thread s. In this formula, which is a generalization of equation (21) of iDymova fc Ruderman 



( 120051 ) . Cfc is analogous to the kink speed of a uniform magnetic tube and is given by cl(z, t) 



2B 2 / fi[p(z,t) + p e ], where p(z,t) must be substituted from equation ([I]). 



3.1. Oscillations of steady threads 



At first let us neglect the observed flows. The norm al modes of a non-flowing promin ence 
thread in the linear, low-/? li mit have been studied by iDfaz. Oliver. &: Ballesterl (120021 ) and 
Dymova &: Rudermanl (120051 ). Such structure supports Alfven and fast waves. The first ones 
cause azimuthal (i.e. torsional) motions and so do not disturb the cylinder axis. Therefore, 
Alfven modes are not the origin of the detected oscillations. Regarding fast waves, the only 
solution that gives rise to a displacement of the cylinder axis is the kink mode, hence it 
provides the most natural interpretation for the observed oscillations. 



Of the two studies just mentioned, that of IDymova fc Rudermanl (120051 ) provides with 
the simplest dispersion relation because of the thin tube approximation. The frequency of 
the kink mode is given by the smallest positive root of their equation (27), in which e = p e j p c , 
c = p p / p c , I = WJ L, and Q = uL/vac, with uj the frequency and v Ac the coronal Alfven speed. 
The Alfven speed in the thread and the evacuated parts of the tube are va p . = e^^VAr and 



c 1/2 v Ac , 



where v\ 



B 2 /(pp ) has been used. According to IDymova fc Ruderman 



(120051 ) the kink mode frequency is practically constant for 0.6 < e < 2 (see their Figure 4). 
This result thus implies that, if the densities in the corona and in the evacuated part of the 
loop are similar, then e = 1 can be take n with out loss of generality. With this approximation 
equation (27) of IDymova fc Rudermanl (120051 1 reduces to 
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tan [n(l - I)] - \lY^ cot 



m > 1 + c 



0. (3) 



The observations of lOkamoto et aJj (120071 ) provide us with two useful quantities: the 



oscillatory period (from which uo is computed) and the thread length (equal to 2W in our 
model). On the other hand, equation fl3} results in a single value of Q once c and I are 
fixed and so it seems that this seismological problem has a unique solution. One must bear 
in mind, however, that the dimensionless variables c and I actually hide five dimensional 
variables (p c , p p , L, W, and va c ) and thus there are infinite solutions. Still some information 
can be extracted if we proceed as follows: let us start by selecting a thread from table [H and 
let us fix a total length of the magnetic tube, for which we have estimated a minimum value 
of 100,000 km, or L > 50, 000 km. Thus, L and W are known and so I can be determined. 
Next, let us fix the density ratio, c = p v j p c , and let us calculate the smallest positive root, fl, 
of equation ([3]). Now the definition of Q can be used to determine the coronal Alfven speed, 
v Ac — ujL/Q = 2ttL/ (£1P), where P is the period, given in tabled) Finally, the Alfven speed 
in the prominence is obtained, va p = c"^ 2 vac- Hence, this procedure gives a unique pair of 
values vac, va p for each pair of L, c; in other words, if we impose a length of magnetic field 
lines and the density ratio, c, for a given thread, then the coronal and prominence Alfven 
speeds are well determined. To display the results of this analysis we keep L fixed and vary 
the density ratio, c, between a value slightly larger than 1 (extremely tenuous thread) and 
1000 (very dense thread), and then plot function of va c for this range of c. Note 

that, although the coronal and thread densities cannot be computed from their respective 
values of the Alfven speed because the magnetic field strength is unknown, the density of 
prominence plasmas is typically two orders of magnitude larger than that of the corona, 
which corresponds to c of the order of 100. 



The results for two of the six threads investigated by lOkamoto et al.l (120071 ) are plotted 
in Figure [2j the other four threads have also been studied and lead to qualitatively similar 
results to those of this figure. The leftmost point of the curves corresponds to c > 1, that 
is, to a thread density slightly larger than the coronal one (this point is not visible because 
it falls outside the range of the vertical axis). As a consequence, the solution va p > va c 
is obtained. When the density ratio is increased from this minimum value we move along 
each curve towards higher coronal Alfven speeds. The left part of all curves displays a 
rapid variation of va p , but for larger values of c the Alfven speed in the thread stabilizes at 
an approximately constant value, in spite of the unrealistically large coronal Alfven speeds 
attained in these plots. From these results it is clear that a lower limit for the Alfven speed 
in the threads can be established. For a length of magnetic field lines of 100,000 km this 
lower bound is between 200 km s _1 and 350 km s _1 , with an exceptional lower value of 120 
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km s^ 1 for thread #6. These values of va p are consistent either with a strong magnetic field, 
e.g. 50 G, and a large density, of the order of 1 — 8 x 10 11 cm -3 , but also with a weaker 
magnetic field and lower density (10 G and 0.4 — 3 x 10 10 cm -3 , for example). 



3.2. Oscillations of flowing threads 

Next the flowing motion of the threads is included in our model and again p e = p c 
is consid ered for simp licity. Equation (j2J) has been integrated numerically using the code 
PDE2D (ISewelll 120051 ) . The large photospheric inertia in front of coronal perturbations is 
taken into account by imposing the boundary conditions u = at z = ±L. The initial 
conditions for u and du/dt must be carefully chosen in order to excite only the global 
transverse motions. Thi s is accomplished by imposin g the normal mode profile at t = 0, 



given by equation (24) of iDymova fc Rudermanl (120051 ) once the frequency is computed from 



equation fl3]); this condition for u(z,t = 0) is accompanied by du/dt(z,t = 0) = 0. The 
normal mode profile is more easily computed when the density is symmetric about the 
center of the tube, and this is the reason for choosing the thread initially placed in the tube 
center at t = in equation ([!]). The duration of the numerical simulations is chosen so as to 
prevent the thread material from reaching the photosphere, and this typically corresponds 
to 4-6 kink mode periods. A uniform grid of 3000 points in the z-direction is used in all 
simulations, which ensures a minimum of 50 grid points in the thread. 

The numerical code has been run for each thread, whose length, flow speed, and oscil- 
latory period are given in table and for the parameter values (length of magnetic field 
lines, density ratio, and coronal and thread Alfven speeds) corresponding to the symbols in 
Figure [2J We first set v = and so expect to obtain the whole tube oscillating in the kink 
mode if the code works well. A representation of u versus z for different times shows that 
the whole thread oscillates in phase, which is consistent with the kink mode being excited 
by the initial perturbation. In addition, the power spectra of u at various fixed positions 
present a strong power peak at a frequency that perfectly matches that of the kink mode. 
We are thus confident about the performance of our numerical code. 

To assess the importance of the flow motion on the oscillatory period the above simu- 
lations are repeated with vq taken from table HJ Again it is found that the tube oscillates 
bodily in phase as it moves at the flow speed, vq. The period is determined from the position 
of the peak in the power spectrum of the signal at a given point and then compared to the 
period given in table [TJ In § 13.11 the values in this table were used to obtain the results 
of Figure [2, so any difference between the observed periods and those from the simulations 
with v o 7^ comes from the thread flow. Our main conclusion here is that all the numerical 
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simulations yield a period of transverse oscillations that is slightly smaller than that of the 
kink mode. The largest influence of the flow on the period is for short magnetic tubes and, 
for a given tube length, the largest deviations from the case Vq — occur for the largest 
density ratios, although it stops increasing for c of the order of 100. For example, taking 
a total length of the magnetic tube of 100,000 km and p p / p c = 200 this difference ranges 
between 3% and 5% for the six threads. For longer magnetic tubes this relative difference 
becomes even smaller. When these numbers are compared with the error bars of the period 
in table [TJ whose minimum value is 6.8%, it turns out that including the flow in the theo- 
retical calculations leads to variations in the period that are undetectable with the cadence 
used to obtain the data. 

To test the above results, o btained in the thin tube app roximation, the linear, ideal 



MHD equations (eqs. (la)-(ld) in lDymova fc Rudermanll2005l . with the replacement of d/dt 



by d/dt + Vod/dz inside the tube) have been solved. The azimuthal dependence exp(im0), 
with m = 1, has been imposed and so a set of two-dimensional time-dependent equations has 
been investigated. No further approximations are made, i.e. the thin tube approximation 
is not used and the flow is maintained in the equations. The obtained results confirm the 
accuracy of the previous calculations. 



4. Discussion 

In this work we have performed a seismological analysis of some oscillating threads in 
an active region prominence observed with Hinode SOT. The horizontal threads flow along 
a path parallel to the photosphere and undergo simultaneous transverse oscillations that are 
interpreted as a signature of the kink mode of the whole magnetic tube in which the thread 
resides. Our numerical simulations prove that when the thread motion is included the global 
transverse oscillations persist and can be initiated by an external impulse. Therefore, the 
most natural explanation for the se transverse thread oscillatio ns is the excitation of the kink 
mode (see also the discussion in Van Doorsselaere et al]|2008l ) , rather than Al f ven waves or 



motions along a helically twisted tube such as put forward by lOkamoto et al.l (120071 ). 



The available data are insufficient to derive well-constrained values of the physical vari- 
ables, although it has been possible to establish a lower limit for the Alfven speed in each 
of the threads. This lower bound comes from the assumption of the same minimum length 
for all threads, namely 100,000 km, but it could be larger if the actual length of the mag- 
netic tube along which the thread flows is larger. To obtain other plasma parameters from 
the thread Alfven speed one must make some assumptions. For example, if a magnetic field 
strength typical of active region prominences is used (50 G), then large densities (1 — 8 x 10 11 



- 8- 



cm 3 ) typical of active region prominences are obtained. 



We have obtained a dispersion of the minimum Alfven speed between 120 km s _1 and 
350 km s -1 that can have multiple interpretations. On one hand, it is possible that not all 
threads have the same length, as has been assumed to derive these values, and that once the 
actual length is used all these values become the same. But on the other hand, one should 
consider the possibility of a highly inhomogeneous prominence, in which the magnetic field 
and density vary in space. This inevitably leads to an inhomogeneous distribution of the 
Alfven velocity in the prominence, which is not surprising given the complex structuring of 
these objects. 

Another important conclusion of this work is the insensitivity of the period to the flow 
velocity, which is a useless parameter in the present seismological problem. This has allowed 
us to simplify our study by considering a static situation. In addition, the thread thickness 
is another parameter of little importance given the much larger extent of the magnetic tube. 

Some simplifications have been made to facil itate our theoretical investigat i on. F irst of 
all, field line curvature has been neglected since iTerradas. Oliver, fc Ballesterl (120061 ) have 
shown that it has little influence on the transverse vibrations of a flux tube, in the context 
of coronal loop oscillations (see also I Van Doorsselaere et al.l 12004 ) . In addition, the effect 
of the gas pressure in the MHD equations has been discarded since the plasma (3 is much 
smaller than unity in the corona and in active region prominences. Another approximation 
has been to neglect non-linear terms in the MHD equations, which is justified since the 
velocity amplitude of transverse oscillations (whose largest value is 22 km s -1 ; see tabled]) 
is much smaller than the Alfven speed both in the threads and in the corona. 

A final comment can be made about the equilibrium model, that does not include 
gravity and therefore does not take into account the necessity of the magnetic structure 
to be able to support the dense prominence mater i al. T hi s effect has only been i n clude d 
in similar flux tube models by iBallester &: Priestl (119891 ); ISchmitt fc Degenhardtl (119951 ); 
Rempel. Schmitt. &; Glatzell (119991 ) . but the resulting structure is unstable. Hence, there 
is currently no prominence model that incorporates all relevant physics and leads to rea- 
sonable (stable) solutions and so using the simpler configuration studied in this paper is 
justified. 
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Table 1. Summ ary of geometric and w ave properties of vertically oscillating flowing 
threads analyzed by lOkamoto et al.1 (120071 ). 2W is the thread length, v$ its horizontal flow 
velocity, P the oscillatory period, V the oscillatory velocity amplitude, and H the height 

above the photosphere 



Thread 


2W (km) 


vo (km s 1 ) 


P(s) 


V (km s- 1 ) 


H (km) 


1 


3600 


39 


174 ± 25 


16 


18,300 


2 


16,000 


15 


240 ± 30 


15 


12,400 


3 


6700 


39 


230 ± 87 


12 


14,700 


1 


2200 


46 


180 ± 137 


8 


19,000 


5 


3500 


45 


135 ± 21 


9 


14,300 


6 


1700 


25 


250 ± 17 


22 


17,200 




Fig. 1. — Sketch of the magnetic and plasma configuration used to represent a flowing thread 
(shaded volume) in a thin magnetic tube. The two parallel planes at both ends of the cylinder 
represent the photosphere. 
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Fig. 2. — Dependence of the Alfven velocity in the th read (va p ) as a functio n of the coronal 
Alfven velocity (va c ) f° r two of the threads studied by lOkamoto et al.l (I2007I ). In each panel, 
from bottom to top, the curves correspond to a length of magnetic field lines of 100,000 km, 
150,000 km, 200,000 km, and 250,000 km, respectively. Asterisks, diamonds, triangles, and 
squares correspond to density ratios of the thread to the coronal gas p p / p c — 5, 50, 100, 200. 



